Nondispersive two-electron Trojan wave packets 
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We demonstrate the existence of stable non-dispersing two-electron Trojan-like wave packets in 
the helium atom in combined magnetic and circularly polarized microwave fields. These packets 
follow circular orbits and we show that they can also exist in quantum dots. Classically the two 
electrons follow trajectories which resemble orbits discovered by Langmuir and which were used 
in attempts at a Bohr-like quantization of the helium atom. Eigenvalues of a generalized Hessian 
matrix are computed to investigate the classical stability of these states. Diffusion Monte Carlo 
simulations demonstrate the quantum stability of these two-electron wave packets in the helium 
atom and quantum dot helium with an impurity center. 

PACS numbers: 31.50.+W, 32.80.Rm, 42.50.Hz, 95.10.Ce 



The direct manipulation of atoms and ions at the quan- 
tum level is currently a flourishing area of physics Q. 
For example, ion traps combined with laser cooling tech- 
niques have been used to create new states of matter 
including ion liquids, Wigner crystals and Bose-Einstein 
condensates |2, y] • The ability to manipulate the quan- 
tum properties of matter directly, as exemplified by these 
advances, is central to the practical development of nan- 
odevices, e.g., quantum dots and microchip traps [4|. 

Recently it has proved possible to create a trap consist- 
ing of a single atom inside of which the quantum behavior 
of an electron can be manipulated 5]. In these experi- 
ments an electron in an excited lithium atom was 
localized in a classical orbit almost indefinitely, neither 
spreading nor dispersing. This "classical atom" was syn- 
thesized by "tethering" the electron using a microwave 
field to which its motion is phase locked. Potential prac- 
tical applications of the technique include Rydberg tag- 
ging in molecular spectroscopy and the preparation of 
stable antimatter atoms 0. 

Although only recently realized experimentally, the 
existence of coherent, non-dispersive one-electron wave 
packets in Rydberg atoms was predicted about a decade 
ago using essentially classical mechanical arguments 
E3 E| • A peculiar property of this type of wave packet 
is that the electron is localized at an equilibrium which 
corresponds to an energy maximum in the noninertial 
frame. Such equilibria are similar to the well known La- 
grange equilibrium points L4 and L5 in the restricted 
three-body problem of celestial mechanics at which, e.g., 
Jupiter's Trojan asteroids are located 0,^1- Because 
their classical dynamical origins are so similar, these wave 
packets have come to be known as Trojan wave packets 

Naturally it is interesting to ask if similar wave pack- 
ets but involving more than a single electron can be syn- 
thesized. In fact, considerable effort has already been 
directed towards this goal. For example, two-electron 
wave packets in barium atoms [13 have been prepared 
and provide a powerful tool in the study of electron cor- 



relation in atoms However, while these wave packets 
are prepared initially in localized radial states using laser 
pulses they eventually disperse, e.g., through a collision 
near the nucleus. In another approach a particular type 
of dynamical stabilization is used to generate nondispers- 
ing wave packets in two-electron atoms. In this case it 
is necessary that one of the electrons is well localized 
relatively close to the nucleus [l^ . 

Here we report the existence of coherent, two-electron, 
non-dispersing wave packets in helium-like atoms (or 
quantum dot helium) which are true analogs of the one- 
electron Trojan wave packets; both electrons follow clas- 
sical orbits in direct analogy with the states prepared 
experimentally in Refs. 0,0. Unlike the one-electron 
Trojan problem, this system is a genuine quantum three- 
body problem and thus represents a more direct analogy 
with the classical restricted three-body problem. Stable 
two-electron equilibria are produced through the simulta- 
neous application of combined circularly polarized (CP.) 
electromagnetic and magnetic fields to the helium atom 
or to quantum dot helium with an impurity center [l^ . 
The equilibria so produced are stable over broad ranges 
of field parameters and the wave packets are, therefore 
true non-dispersing coherent-like states. 

We demonstrate that, in the two-dimensional limit, 
e.g., as in quantum dot helium, it is possible to create 
a variety of stable two-electron, non-dispersive, Trojan 
wave packets. However, there exists only a single stable 
three-dimensional configuration in the helium atom it- 
self. These latter wave packets are actually two-electron 
examples of the electronic coherent states sought after 
by Schrodinger in the hydrogen atom. In particular 
these wave packets follow the "double-circle" orbits orig- 
inally discovered by Langmuir in the helium atom. 
We also present diffusion Monte Carlo simulations which 
confirm directly the quantum stability of these states. 
Further, we show that it should be possible to produce 
certain of these states directly in quantum dots using 
field strengths which are currently accessible in the lab- 
oratory. 
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where qk = x,y,z\k = 1,2,3. Setting the velocities to 
zero produces the ZVS which, in the present case, has 
the form 
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Vzvs = VzvS! + Vzvs 2 + V r 12 



(4) 



FIG. 1: Possible rotating configurations of two electrons in 



a helium atom or quantum-dot helium in combined magnetic 
and CP fields. Type I: Langmuir configuration, type II: trans- 
verse configuration, type III: collinear configurations (a and 

b). 



In a coordinate system rotating with the CP. field and 
assuming an infinite nuclear mass the Hamiltonian for 
the helium atom interacting with a CP. field and a mag- 
netic field perpendicular to the plane of polarization is, 
in atomic units, 



where, 
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H = H 1 +H 2 + — 
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withz = 1,2; n = V x i +Vi + z i and P* = (PxuPynPzi) 
are the coordinate and momentum vectors of each elec- 
tron; uo and £ are the scaled CP. field frequency and 
scaled field strength, respectively [Tjj ; the cyclotron fre- 
quency in scaled units is ±1. With a re-interpretation 
and re-scaling of parameters this is also the Hamiltonian 
of quantum dot helium with an off- axial impurity cen- 
tre in a magnetic field Equilibria of the classical 
Hamiltonian can most easil y b e found by constructing a 
zero- velocity surface (ZVS) [lfl|, which shares some prop- 
erties with a potential energy surface (PES) but contains 
addition terms due to the centrifugal and Lorentz forces 
[HI EI . An extended discussion of the construction of 
the ZVS in atomic systems is given in - briefly the 
ZVS is obtained by re- writing the Hamiltonian in terms 
of velocities rather than momenta. This gives 
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The ZVS reduces to a true PES when uj = 1/2 , i.e., 
when the paramagnetic term in eq. (2) vanishes. As 
with a PES equilibria of the motion are then obtained 
as extrema of the ZVS (maxima, minima and saddles) 

nam. 

We find three possible types of equilibrium point whose 
stability depends on the particular values of the parame- 
ters. Type I: This configuration, illustrated in Fig. 1 (I), 
corresponds to the unstable "two-circle" orbits discovered 
by Langmuir |l6| in which the (classical) electrons occupy 
parallel orbits located above and below the plane of the 
nucleus with longitudes of 30°. Stabilization is possible 
only through the application of external fields. Type II: 
Illustrated in Fig. 1 (II), this configuration has the two 
electrons orbiting in the same orbit in a plane contain- 
ing the nucleus. The angle subtended by the electrons 
at the nucleus can be arbitrary. This configuration can 
only be stable in two spatial dimensions e.g., a quantum 
dot. Type III: This configuration, of which two variants 
are possible - Ilia and Illb in Fig. 1 - is unstable in the 
helium atom itself but stable in 2-dimensional quantum 
dot helium where it appears as variant Ilia. Both elec- 
trons lie in the plane of polarization and, remarkably, lie 
on the same side of the nucleus. 

The equilibrium corresponding to the Langmuir con- 
figuration has the geometry of an equilateral triangle; 
i.e., the three particles form an equilateral triangle whose 
sides are of length a (see Fig. 1 a and Fig. 2). The rela- 
tionship between the field parameters and the size of this 
triangle is defined by the following cubic in a; 
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Note that the presence of external fields does not change 
the geometry of the Langmuir configuration. Figure 2 
shows the dependence of a on the electric field strength 
for selected frequencies uj. We perform an extended sta- 
bility analysis of the system around these trajectories, i.e. 
we study small oscillations |20(, in the rotating frame. 
The generalized Hessian matrix (stability or monodromy 
matrix) for a system of n-particles in the rotating frame 
can be written as 



where 
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FIG. 2: Dependence of the length of the sides of the equi- 
lateral triangle described by the three charges in the rotating 
frame on the strength of the CP field for various rotation fre- 
quencies. Note that for uo c = — 1 the system is bistable anc 
only the larger solution corresponds to the stable trajectory 
When cdc = — 1 there is no critical field and only a singk 
solution E3. 



where sgn(ai) = l,a^ = p^ sgn(ai) = — l,a^ = X{ anc 
H e is the phase space defined Hessian matrix 



(H e ) ai 6i 



d 2 H n 
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dctidbi 

is the harmonic Hamil- 



<H = Qi,Pi, h = Qi,Pi and H ns 
tonian around the orbit 0, |2l| . 

The generalized Hessian matrix can be written explic- 
itly as 
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Note that the right upper block is the normal Hessian 
matrix of the stationary mechanical system dj [2£| and 
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Figure 3 shows the stability island as a function of scaled 
coordinates |S| for the cases uo c = 1 (lower half) and 
ou c = —1. Remarkably, unlike for Trojan wavepackets 
both the circularly polarized field and the magnetic field 
are necessary to stabilize the Langmuir trajectories. One 
may notice an interesting v-shaped valley which is similar 
to the "negative mass" case of Ref. . 

The quantum wave function for the monodromy matrix 
for the stable trajectory [2]] may be written 



'Y2 i 7 AijXiXj 



(12) 



where An are always real and Aij, i j imaginary to 
guarantee normalizability of the quantum wavefunction 
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FIG. Stability diagram (scalecl units) of the disturbed 
Langmuir trajectory for ant i- centrifugal Lorentz force (upper 
half, uj c — —1) and co- centrifugal (lower half, uj c = +1) as 
a function of scaled rotation frequency and the electric field. 
The black regions corresponds to stable trajectories. Note 
that no magnetic field cj c < can lead to stability in case —1 
and magnetic field weaker than rotation co c < co is permitted 
to obtain stability in case +1. This suggests a maximum of 
the ZVS and Paul trap-like (Trojan) stabilization. 



in the rotating frame. Because the wavefunction is local- 
ized around the stable point in the configuration plane 
the contraction (integration) of diagonal of the density 
matrix (x)ijj(x) over a single electron variable leads to 
the single electron density localized around electron equi- 
libria and corresponds to a nondispersing two-electron 
wavepacket moving in the laboratory frame. The detailed 
calculations are quite cumbersome and will be given else- 
where. As one can see from Fig. 3 the system never sta- 
bilizes for magnetic fields weaker than the electric field 
frequency for the case of anti-centrifugal Lorentz force 
and may stabilize for such fields for faster rotations for 
the co-centrifugal Lorentz force. 

In order to check for the existence of the quan- 
tum states themselves we solved the time dependent 
Schrodinger equation using the Diffusion Monte Carlo 
method (23, l24j |. Fig. 4 shows density plots of the 
wave function of each electron. The probability den- 
sity is clearly localized around the equilibrium points. 
The wavepackets are actually approximate eigenstates of 
the Hamiltonian in the rotating frame and, therefore, do 
not disperse in the inertial frame. These eigenstates are 
localized around a classical equilibrium point whose lo- 
cal spectrum is an almost-harmonic ladder of coherent 
states. That is, the eigenstate of the atom in the fields 
itself behaves as a wavepacket evolving according to the 
classical equations of motion. For the large field values 
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FIG. 4: L Hftff&8it>ftO2Monteo0aPlo5 simulations of Langmuir 
wavepackets in the rotating frame for Q c = 0.0370 a.u., 
£ — 0. 1235a. u., and Q = Q c /2. The density of points re- 
flects the wavefunction itself not probability density [lgL |23|. 
The wavepackets are approximate eigenstates of the Hamilto- 
nian in the rotating frame and follow a circular orbit in the 
inert ial frame. 

employed here the Langmuir wavepackets actually cor- 
respond to the ground state of eq. (2). This is similar 
to the case of the original Trojan wavepackets 7] which 
are the ground states of the locally harmonic Hamilto- 
nian obtained by expansion around the Trojan equilibria 
which are, in that case, energy maxima |8|. We neglect 
finite nucleus mass effects here but on the timescales of 
interest, as in the case of one-electron Trojan wavepack- 
ets, these effects will be negligible 25]. 

In conclusion we have demonstrated that Langmuir 
trajectories in magnetic and circularly polarized fields are 
stabilizing for certain parameter regions when the mag- 
netic field is parallel to the rotation axis and the CP. 
field is perpendicular. Corresponding quantum states ex- 
ist. This is also true for two-dimensional quantum dot 
helium. These regions are not possible to predict from 
purely analytical considerations and require extensive nu- 
merical searches in parameter space. For example, for a 
quantum dot of radius 100 nm (Type Ilia configuration) 
the equilibria are at 62.72 nm and 98.00 nm from the 
nucleus. To achieve this configuration a magnetic field of 
5 T is applied and the impurity with the effective charge 
Zeff =0.008 e displaced from the center of the parabolic 
dot by 98 nm generates effective "microwave" field of fre- 
quency 548 GHz and strength 4.593 kV/m. 
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